Is the Casimir effect relevant to sonoluminescence? 
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The Casimir energy of a solid ball (or cavity in an infinite medium) is calculated by a direct 
frequency summation using the contour integration. The dispersion is taken into account, and the 
divergences are removed by making use of the zeta function technique. The Casimir energy of a 
dielectric ball (or cavity) turns out to be positive, it being increased when the radius of the ball 
decreases. The latter eliminates completely the possibility of explaining, via the Casimir effect, the 
sonoluminescence for bubbles in a liquid. Besides, the Casimir energy of the air bubbles in water 
proves to be immensely smaller than the amount of the energy emitted in a sonoluminescent flash. 
The dispersive effect is shown to be inessential for the final result. 
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1. Sonoluminescence being observed during more than half century [[lj has not received satisfactory explanation 
yet. As known this phenomenon represents the emission of visual light by spherical bubbles of air or other gas injected 
in water and subjected to an intense acoustic wave in such a way that the radius of bubbles changes periodically. In 
the last years of his life Schwinger proposed that the bases of sonoluminescence is formed by the Casimir effect. 
While changing the size of bubbles the zero point energy of the vacuum electromagnetic field (the Casimir energy) 
of a cavity in a dielectric medium changes too. According to Schwinger, it is these changes of the electromagnetic 
energy that are emitted as a visual light in sonoluminescent flashes. In Schwinger's calculations the Casimir energy 
for the configuration in hand proves to be of the same order as the energy of the photons in an individual flash 
(~ 10 MeV). Other authors obtained results both consistent with Schwinger's calculation JjJ and differing from it 
in 10 orders [Qj5). This disagreement is basically due to different methods used for removing the divergences in the 
problem under consideration. 

In the present note the calculation of the Casimir energy of a dielectric ball placed in an endless dielectric medium 
(or cavity in this medium) is carried out under following conditions. In the first place a realistic description of 
dielectric properties of media is used which takes into account dispersion |(| . On the other hand the most simple and 
reliable method for removing the divergences, the zeta function technique, is applied. Till now these conditions were 
not combined in studies of the problem in question. 

2. When calculating the Casimir energy we shall use the mode- by- mode summation of the eigenfrequencies of the 
vacuum electromagnetic oscillations by applying the contour integration in a complex frequency plane ■ Consider 
a ball material of which is characterized by permittivity E\ and permeability The ball is assumed to be placed in an 
infinite medium with permittivity e 2 and permeability \x 2 . For this configuration the frequencies of transverse-electric 
(TE) and transverse-magnetic (TM) modes are determined by the equations js] 

A™(aw) = V £ iM2 s' l (kia)ei(k 2 a) - y/e^ siikx^e^a) = 0, (1) 
A™ (aw) = V £ 2Mi s' l (k 1 a)ei(k 2 a) - y/e^ s^k^e'^a) = 0, (2) 

where §i(x) = \/ttxJ2 J l+l / 2 {x) and ei(x) — y^nx/2 ; 2 (x) are the Riccati-Bessel functions, ki — yj£i\ii u>, i = 1,2 
are the wave numbers inside and outside the ball, respectively; prime stands for the differentiation with respect to 
the argument (k\a or k 2 a) of the Riccati-Bessel functions. 
As usual we define the Casimir energy by the formula 

£ = 5X>.-S.), (3) 

s 

where u> s are the roots of Eqs. (|l|) and (||) and ui s are the same roots under condition a — > 00. Here s is a collective 
index that stands for a complete set of indices specifying the roots of Eqs. (jl|) and (||): s = {I, m, n} I = 1, 2, . . . ; m = 
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— (I + 1), —I, . . . , I + 1, n = 1,2, The roots of Eqs. (|j) and (^) do not depend on the azimuthal quantum number 

to. Therefore the corresponding sum gives a multiplier (21 + 1). Further we use the principle of argument theorem 
from the complex analysis in order to present the sum over n in terms of the contour integral. As a result Eq. (||) can 
be rewritten as follows: 



i=i 



E, = 
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where the contour C surrounds, counterclockwise, the roots of the frequency equations (|l]) and (Eh in the right half- 
plane. This contour can be deformed into a segment (— iA, iA) of the imaginary axis and a semicircle of radius A with 
A — ► oo. In this limit the contribution of the semicircle into the integral (|j) vanishes with the result || 
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(5) 



x [y/e^m^ ((s'i(qi) e i(l2)) 2 + (sz(gi)ej(q 2 )) 2 ) - (eiM2 + ^Mi^Kft^fejeK^eJO^)] 

where % = -JsifM y, i — 1, 2 and s;(z), e;(z) are the modified Riccati-Bessel functions: s/(z) = (ttz /2) 1 / 2 I ll (z), e;(z) = 
(2z/ir)V 2 K v (z), i/ = 1 + 1/2. 

Further we will content ourselves by examining the case when both the media are nonmagnetic fii = fi2 = 1 and 
their permittivities £\, £2 differ slightly. In view of this we can put in Eq. (^J) q\ = 02 keeping in remain e\ and £2 
exactly. It gives 



E, 



1 + 1/2 



dy\n{l-e [(si(y)ei(y))f} , 



(6) 



v^T + v^2 

Now we are going to take into account the effect of dispersion considering the parameter £ 2 in Eq. (||) as a function 
of y = auj/i. Justification of the mode- by-mode summation method in applying to dispersive and absorptive media 
has been considered in M. For definiteness we put Ex = 1 + 5, £2 = 1, S <C 1, then £ 2 ~ S 2 /16. We substitute 5 by 



5{y)=S /[l + {y/vy f], v = 1 + 1/2, 



(7) 



where Sq is a static value of S(y) and the parameter yo is determined by a "plasma" frequency loq: yo = aujQ. The 
function describing dispersion in Eq. (Q) is a standard one [one-absorption-frequency Scllmcir dispersion relation] 
except for its dependence on We have introduced this dependence in order to be able to use the zeta function 
technique below. This complication does not contradict the main goal pursued by using this function, namely, it 
should simulate crudely the behaviour of 5(y) at large y. As known |l(J, the general theoretical principles lead to the 
following properties of the function e(u>) in the upper half-plane to. On the imaginary axis u> = iy, y > the function 
e(iy) acquires real values, and with increasing y it steadily decreases from the static value 1 + So > (for dielectrics) 
to 1. Obviously formula (R) meets these requirements. 

Substituting (0) into (^J) and making use of the uniform asymptotic expansion for the modified Bessel functions |Tl| ] 
when I — ► 00 one obtains 



Ei ~ -JL(*°) /l(acJo ) 
i->oo 64a V 4 ' J V ' 
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We carry out the summation of the partial Casimir energies 
taking into account asymptotics (eh 



* = E^ = E 
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with the help of the zeta function technique JlJ] 
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7 = 1 



/i(awo)[C(0,l/2) 



Here — Ei + (3/64a) (So/ 4) /i(cwo) is the renormalized partial Casimir energy, f(s, g) is the Hurwitz zeta function. 



As C(0, 1/2) = 0, we get for the Casimir energy (|12| 
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With allowance for (0) one can obtain the estimation for the sum X)i=i &i 
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Thus the Casimir energy of a dielectric ball is 
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(13) 



(14) 



(15) 



dispersion resulting only in the positive functions /,(aa;o), i = 1, 2, 3. When z increases the functions fi(z) approach 1 
(see Fig. 1), and (15) turns into the expression for the partial Casimir energy of a solid ball without dispersion ||. 

Considering the behaviour of the functions fi(z) (see Fig. 1) one concludes that the main contribution, with a few 
percents accuracy, gives the first term in braces in Eq. (nq) with the result 



E 



3 
64^ 



/i(aw ). 



(16) 



Obviously the change of the energy sign or a considerable increasing its magnitude due to the dispersion effect [[L3[ is 
out of the question. 

Let us estimate the value of fi(awo). The parameter ujo can be determined by demanding that at this frequency 
the photons do not 'feel' the interface between two media. This condition will be certainly met when the wave length 
of photon is less than the interatomic distance in media d ~ 10~ 8 cm. Actually it is the condition of applicability of 
the macroscopic description of dielectric media fiofl . Sonoluminescence is observed with the air bubbles in water [[l], 
the radius of bubbles being a ~ 10 -4 cm. Hence it follows that aojo — a/d = 10 4 and /i(10 4 ) = 0.999 . . .. Thus the 
allowance for the dispersion in calculating the Casimir energy of a dielectric ball (or spherical cavity in a slab of a 
dielectric) practically has no effect on the final result. 
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Certainly the real picture of dispersion in the whole frequency range < ui < oo for any dielectric, including water, 
is exceedingly complicated and cannot be described by a simple equation ^ with a single parameter u>q. As known 
absorption of the electromagnetic waves in water and, as a consequence, their dispersion take place already in the 
radio frequency band. Putting in this case A ~ 10 4 cm, we obtain awo ~ 1 and /i(l) = 0.729 . . .. From here one can 
infer that the effective value of own should be less than 10 4 . In order for a more precise evaluation of this parameter 
to be done a more detailed consideration of the dispersion mechanism is needed. Obviously this may lead only to 
diminution of the absolute value of the Casimir energy. However this issue is beyond the scope of the present paper 
for the main conclusion (see below) does not depend on this point. 

It is worth noting two peculiarities of the final formula ( |l6| ) . When the radius of the bubble decreases its Casimir 
energy increases. This behaviour is completely opposite to one needed for explanation of sonoluminescence (as known, 
emission of light takes place at the end of collapsing the bubbles in liquid). Besides, this energy is immensely smaller 
than the amount of energy emitted in a separate sonoluminescent flash (~ 10 MeV). Actually taking a = 10 -4 cm 
and 5q = 3/4 (water) we arrive at E ^ 5 • 10 -3 eV. 

Thus the results of this paper unambiguously testify that the Casimir effect is irrelevant to sonoluminescence. 

This work was accomplished with financial support of Russian Foundation of Fundamental Research (Grant 97- 
01-00745). 
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1'ig.l. Functions fi(z), i 1 ,2,3 defining the expansion (8). 

When z increases the functions approach their limiting 
value equal to 1. 



